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Using our new Post-Newtonian SPH (smoothed particle hydrodynamics) code, we study the final 
coalescence and merging of neutron star (NS) binaries. We vary the stiffness of the equation of 
state (EOS) as well as the initial binary mass ratio and stellar spins. Results are compared to those 
of Newtonian calculations, with and without the inclusion of the gravitational radiation reaction. 
We find a much steeper decrease in the gravity wave peak strain and luminosity with decreasing 
mass ratio than would be predicted by simple point-mass formulae. For NS with softer EOS (which 
we model as simple T = 2 polytropes) we find a stronger gravity wave emission, with a different 
morphology than for stiffer EOS (modeled as P = 3 polytropes as in our previous work). We also 
calculate the coalescence of NS binaries with an irrotational initial condition, and find that the 
gravity wave signal is relatively suppressed compared to the synchronized case, but shows a very 
, significant second peak of emission. Mass shedding is also greatly reduced, and occurs via a different 

mechanism than in the synchronized case. We discuss the implications of our results for gravity wave 
astronomy with laser interferometers such as L1GO, and for theoretical models of gamma-ray bursts 
, (GRBs) based on NS mergers. 
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I. INTRODUCTION AND MOTIVATION 
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Coalescing neutron star (NS) binaries are among the most important sources of gravitational radiation for detection 
by LIGO 0, VIRGO [§, GEO |,||, and TAMA g. These interferometers will be most sensitive to gravity waves 
in the frequency range from ~ 10 Hz to ~ 300 Hz, corresponding to the last few thousands of orbits prior to final 
merging. During final merging, the characteristic frequency is ^ 1 kHz, requiring the use of special narrow-band 
detectors || , which are currently being tested on the GEO 600 interferometer in Germany || . The detection of these 
^ merger wave forms, combined with theoretical knowledge about the hydrodynamics of the merger process, could yield 
I \ detailed information about the behavior of the NS fluid, and, in particular, the EOS of nuclear matter at high density. 
J-Hi This is the second of a series of papers in which we attempt to develop this theoretical knowledge using 3D numerical 
ZV\ hydrodynamics calculations of NS mergers in post-Newtonian (PN) gravity. 

The first Newtonian calculations of binary NS coalescence were performed more than 10 years ago by Nakamura 
and Oohara using a grid-based Eulerian finite-difference code Rasio and Shapiro (Ref. ||], hereafter RS1-3 or 

collectively RS) later used Lagrangian SPH (smoothed particle hydrodynamics) calculations to study the stability 
properties of close binary NS configurations in strict hydrostatic equilibrium, and follow the evolution of unstable 
systems to complete coalescence. RS demonstrated for the first time that, even in Newtonian gravity, an innermost 
stable circular orbit (ISCO) is imposed by global hydrodynamic instabilities, which can drive a close binary system 
to rapid coalescence once the tidal interaction between the two NS becomes sufficiently strong. Since then many 
other groups have performed Newtonian calculations of binary NS coalescence, concentrating on various aspects of 
the problem. Most recent calculations have used either SPH (see, e.g., Jl^ Q)' or PPM (the Eulerian, grid-based, 
piecewise parabolic method; see, e.g., |lq-|l9|]). Refinements have included the addition of terms to approximate 



energy and angular momentum losses to gravitational radiation (see [|^|P1|17 ) and more detailed treatment of the 
nuclear physics (e.g. fli~3| , |l8| ], with applications to models of gamma-ray bursts). 

The first calculations to include the lowest-order post-Newtonian (1PN) corrections to Newtonian gravity, as well 
as the lowest-order dissipative effects of the gravitational radiation reaction (2.5PN) were performed by Oohara and 
Nakamura using an Eulerian grid-based method (^0|. More recently, Faber and Rasio ( &J, hereafter Paper 1), and 
Ayal et al. p2[ , have performed SPH calculations that included all 1PN and 2.5PN corrections, using a PN formalism 
developed by Blanchet, Damour, and Schafer ( |23|, hereafter BDS). 

These calculations revealed that the addition of 1PN corrections to the hydrodynamic equations can have a marked 
effect on the gravity wave signal. In Paper 1, we presented a detailed comparison between two coalescence calculations 
(for initially synchronized binaries containing two identical stars modeled as simple r = 3 polytropes), one including 
both 1PN and 2.5PN corrections, and the other including only Newtonian gravity plus radiation reaction (2.5PN 
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corrections). We found that the 1PN effects produce several luminosity peaks in the gravity wave emission, whereas 
the Newtonian merger shows a simpler signal with a single peak followed by exponential damping of the emission. 
In addition, 1PN effects accelerate the final coalescence and produce large tidal lag angles for both NS just prior 
to contact, leading to a more "off-center" collision. This produces a complex pattern of quadrupole oscillation and 
differential rotation in the merger remnant, accompanied by a characteristic modulation of the gravity wave signal. 

Several groups have been working on fully general relativistic (GR) calculations of NS mergers, combining the 
techniques of numerical relativity and numerical hydrodynamics in 3D pi[-|2q]. However, this work is still in the 
early stages of development, and only preliminary results have been reported so far. Obtaining accurate gravitational 
radiation wave forms from full GR simulations is particularly difficult, since the waves must be extracted at the outer 
boundaries of large 3D grids extending out into the true wave zone of the problem. However, fully GR calculations are 
essential for addressing the question of the stability to gravitational collapse and ultimate fate of a merger. Using a 
GR formalism developed by Shibata (2^] , Shibata and Uryu j2^] find that collapse to a black hole on a dynamical time 
scale during final coalescence does not happen for NS with stiff EOS and realistic parameters. Only with unrealistically 
compact initial NS models (starting already very close to the maximum stable mass) does collapse to a black hole 
occur on the merger time scale. This confirms our expectations from the discussion in Paper I: for realistic NS with 
stiff EOS, PN calculations of binary coalescence can provide results that remain at least qualitatively correct all the 
way to a complete merger. 

In this paper, we continue the work begun in Paper 1, using our PN SPH code to study the coalescence of NS 
binaries. Here, we relax some of the simplifying assumptions made in Paper I in constructing the initial models 
for NS binary systems. We no longer assume the binary system to contain two identical stars, and perform several 
calculations for binary mass ratios q ^ 1. We also allow for a nonsynchronized initial configuration, instead of the 
rigidly rotating binary system assumed in Paper 1. Finally, we consider NS with somewhat softer EOS, but still 
modeled as simple polytropes. In the rest of this section, we discuss the theoretical and observational constraints on 
these various assumptions. 

Mass measurements for NS derived from the timing of relativistic binary radio pulsars are all consistent with 
a remarkably narrow underlying Gaussian distribution with Mns = 1.35 ± 0.04 Mq f27|) . The largest observed 
departure from q = 1 in any known binary pulsar with likely NS companion is currently q — 1.386/1.442 = 0.96 
for the Hulse-Taylor pulsar PSR B1913+16 [§§. For comparison, q = 1.349/1.363 = 0.99 for PSR B2127+11C 
HI, and q = 1.339/1.339 = 1 for PSR B1534+12 f^Q. Although the equal-mass case is clearly important, one 
should not conclude from these observations that it is unnecessary to consider coalescing NS binaries with unequal- 
mass components. Indeed, it cannot be excluded that other binary NS systems (that may not be observable as 
binary pulsars) could contain stars with significantly different masses, especially considering the different histories 
of the primary and secondary in most binary NS evolution scenarios. Theoretically, the predicted NS masses from 
stellar evolution and supernova calculations span a considerably wider range than observed Moreover, Newtonian 
calculations of binary NS coalescence (RS2, |ll]]) have shown that even very small departures from q = 1 can drastically 
affect the hydrodynamics. A small asymmetry in the initial system grows rapidly during the merger, and there is 
almost no disruption of the primary (more massive star), which simply sinks to the center of the merger remnant 
(RS2). A more extreme asymmetry, such as q ~ 0.5, can lead to a brief episode of mass transfer from the secondary 
(less massive star) onto the primary, followed by a sudden widening of the binary system, rather than a merger (RS2). 
The same behavior is seen in hydrodynamic simulations of coalescing binaries containing a NS in orbit around a more 
massive black hole (see H^Q). In this paper, we will not consider such extreme mass ratios, but will focus instead 
on covering well a narrower but realistic range of values from q = 0.8 to q — 1 (For comparison, previous, purely 
Newtonian calculations have been performed for q = 0.5 and q = 0.85; see RS2, |ll|). 

Because of the low viscosity of the NS fluid, the tidal synchronization time for NS binaries is always considerably 
longer than the orbital decay time scale. Therefore, it is expected that the NS should be spinning slowly prior to 
final coalescence ]34| . In the limit of nonspinning stars at large separation, the binary configuration just prior to final 
coalescence will be irrotational, i.e., the fluid motion inside each star will have zero vorticity, although it can have a 
small amount of angular momentum ( |^|-|3(|). In contrast, in many previous studies of NS binary coalescence (and 
in Paper 1), synchronized initial conditions have been used for simplicity (RS2, |^ J| |l^jL6| , pO[ ) . For a synchronized 
binary, the entire mass of fluid is in rigid rotation, and an equilibrium initial configuration can be constructed easily 
using relaxation techniques in the corotating frame of the binary (where the problem reduces to hydrostatics). These 
relaxed initial conditions, while somewhat unphysical, are much less susceptible to spurious oscillations of the fluid 
around the correct quasi-equilibrium configuration (RS). In this paper, we will present a method for constructing 
accurate irrotational initial conditions, and we will compare the results of PN calculations for the coalescence of 
irrotational NS binaries to those obtained in Paper 1 for initially synchronized NS binaries. 

Here, as in Paper 1, we model NS with a stiff EOS as simple polytropes, i.e., the pressure P is a function of the 
density only, of the form P = kr^ , where k is called the polytropic constant, T is the adiabatic exponent, and r» is the 
rest mass density of the fluid (we adopt the notation of BDS), not to be confused with the total mass-energy density 
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of the fluid §7|. 

A polytropic EOS provides a convenient approximation to many different realistic NS EOS calculated under various 
assumptions about the still uncertain microphysics above nuclear density. Higher values of T represent stiffer EOS, 
while lower values represent softer EOS. In addition, it allows for direct comparison to previous Newtonian work, 
which was done mainly for polytropes. The most commonly adopted polytropic EOS for NS uses T = 2 (RSI, 
|^0| , ^5|Jl6| ) . Newtonian polytropes with r = 2 have a radius that is independent of mass, a well-known property 
satisfied approximately by all realistic NS models ( |p8| ). Several tabulated NS EOS are found to have effective 
adiabatic exponents that are close to this value, or slightly higher. In particular, the Lattimer-Swesty EOS used 
in several previous calculations |lli| , |l4| , |l7j -|l9[ has an effective adiabatic index T ~ 2.3 for a 1.4M NS, assuming the 
lowest available value of the nuclear compressibility (K = 180 MeV). The latest microscopic NS EOS, constrained 
by nucleon scattering data and the binding of light nuclei, and incorporating three-body forces, are even stiffer, with 
r ~ 3 (see, e.g., Q for a summary, Q for a recent update, and |3(| for polytropic fits). In this paper, we will present 
and compare calculations performed for r = 3 (as in Paper 1), and for T = 2. We feel that this range corresponds 
reasonably to the theoretical uncertainty in recent determinations of the NS EOS. We will not consider in this paper 
the much softer EOS that may result from even more uncertain considerations of exotic states of matter such as pion 
condensates or strange quark matter J42|| . Note that PN expansions are more difficult to do consistently with softer 
EOS, since the stars are more centrally condensed and thus the central gravitational field is much stronger. There 
is some advantage to using softer EOS within the framework of the BDS formalism, however, since smaller values 
of r help keep several 1PN quantities small (see Paper 1, Eqs. A10-12). A key difference between stiffer and softer 
EOS is captured by our range of adiabatic exponents. Indeed, previous Newtonian work has shown a significant 
difference between T = 2 and T = 3, since polytropes with T > 2.4 can support a stable, rapidly rotating triaxial 
configuration, which would lead to a persistent gravity wave signal, whereas those with T < 2.4 generally damp out 
to oblate spheroids, which have no time-dependent quadrupole moment and therefore cannot radiate gravity waves 
(RS2). The existence of a long-lasting exponentially damped gravity- wave emission tail would have a significant effect 
on the integrated power spectrum of the gravitational radiation signal from binary coalescence. 

The outline of our paper is as follows. Section II presents our numerical methods, including a brief description of our 
PN SPH code, and the procedure used to construct irrotational initial conditions (The construction of synchronized 
initial conditions and the SPH method are discussed in more details in Paper 1). Additionally, we describe the 
parameters and assumptions for all runs quoted in this paper. Section III presents our numerical results, focusing 
in turn on the effects of varying the EOS, the mass ratio, and the initial spins, as well as the differences between 
Newtonian runs and PN runs. Motivation for future work as well as a brief summary of our results to date are 
presented in Section IV. 



II. NUMERICAL METHODS AND SUMMARY OF CALCULATIONS 



A. The PN SPH code 



All our calculations were performed using the post-Newtonian (PN) smoothed particle hydrodynamics (SPH) code 
described in detail in Paper 1. It is a Lagrangian, particle-based code, which has been adapted to the PN formalism 
of BDS. As in Paper 1, we use a hybrid 1PN/2.5PN formalism, in which 1PN corrections are artificially scaled down, 
so as to remain numerically tractable, while 2.5PN effects are included for realistic NS parameters. Throughout this 
paper we use units such that G = M = R = 1, where M and R are the rest mass and radius of the primary NS. In 
these units, there are then two different values of the speed of light c corresponding to the 1PN and 2.5PN corrections. 
As in Paper 1, we take C2.5PN = 2.5, which gives GM / Rc\ bPN = 0.16, corresponding to realistic NS parameters, but 
we take cipat = 4.47, giving GM / Rc\ PN — 0.05, to keep all 1PN quantities sufficiently small. 

Models for the primary NS with a T — 3 polytropic EOS (i.e., P = kr^) were taken from Paper 1. To construct 
models for the secondary NS when q ^ 1, we kept the same value of the polytropic constant k determined for a 
primary with GM/Rc 2 — 0.05, but we computed models with a lower mass using a relaxation technique. Note that, 
in contrast to their Newtonian counterparts, PN polytropes do not obey a simple power-law mass-radius relation, and 
therefore models for the secondary cannot be obtained simply by rescaling models for the primary. Additionally, we 
computed a series of T = 2 polytropic models for single NS with values of GM/Rc 2 between 0.01 and 0.05, as well as 
lower-mass models based on the primary model with GM/Rc 2 = 0.05. While Newtonian polytropes with T = 2 have 
a radius independent of the mass (for fixed k), the radius of a PN r = 2 polytrope increases with decreasing mass. 
The properties of the T = 2 PN polytropic sequence are summarized in Table | and the properties of NS models with 
lower masses are listed in Table [|[ 
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To set up synchronized initial conditions, we used the same relaxation scheme described in Paper 1 to create 
accurate equilibrium configurations. The two NS are initially placed along the x axis, with a binary separation tq 
(distance between the centers of mass of the two stars), and the binary center of mass at the origin. The orbital 
motion is assumed to be in the x — y plane, and the orbital angular momentum vector in the positive z direction. 
We assume that all SPH particle velocities Vj vanish in the corotating frame of the binary, and we use an iterative 
method to find the proper force on each particle in equilibrium, holding the binary separation fixed. As discussed in 
Paper 1, there are subtle differences between the particle velocities Vj and momenta Wj, which must be treated with 
care (see Paper 1, Eqs. A23-24, and Eq. B8). The angular velocity f2 of the corotating system is calculated from the 
centripetal accelerations of the two components of the binary, 



Here v*- 1 ) and v' 2 ' are the center-of-mass velocities of the primary (located on the positive x side) and secondary 
(located on the negative x side), respectively. As in Paper 1, we set the initial velocity field of the fluid to be v x = — fiy, 
v y = Qx when the dynamical run is started (in the inertial frame). 

B. Irrotational initial conditions 

Since there is no simple way to relax an irrotational binary configuration, to strict equilibrium, we used the results 



of Lombardi, Rasio, and Shapiro ( 43 1, hereafter LRS) to construct approximate initial conditions. LRS calculated 
PN equilibrium solutions for irrotational binary NS with polytropic EOS, assumed to have self-similar ellipsoidal 
density profiles, with the density as a function of radius given by the 1PN expansion of the Lane-Emden equation. 
Approximate solutions were determined by minimizing the total energy of the binary configuration, including 1PN 
terms. The resulting NS models are the compressible, PN analogues of the classical Darwin-Riemann ellipsoids for 
incompressible fluids J35],[l4|,[l5| . 

In this paper we will only consider a single irrotational binary model, with q = 1 and r = 3. To construct the 
initial binary model, we adopt the equilibrium shape of the stars determined by LRS, but in combination with the 
internal structure of our single star models (calculated previously for synchronized initial conditions) . Since our PN 
runs were started from a rather large initial separation, ro = 4.0 R, where the tidal deformations are still extremely 
small, we find that the initial oscillations that result from small departures from true equilibrium are negligible. For 
an EOS with T = 3 (n = 0.5, where n = l/(r — 1) is the polytropic index), an initial separation of 4.0 R, and a 
compactness GM/Rc 2 = 0.05, we find from Table III of LRS that the axes of the ellipsoidal figure for each star are 
given by a±/R ~ 1.02, 02/01 — 0.96, and 03/01 ~ 0.96. Here ai is along the binary axis (x-direction) , 02 is in the 
direction of the orbital motion (y-direction) , and 03 is along the rotation axis (^-direction). 

The initial velocities of all SPH particles are calculated from Eqs. 18-20 of LRS, which give 

v x = -Sly(l - -^-2) (2) 
a\ + Oj 

% = - 72^2)- (3) 

ll^ ~r Oo 

It is easy to verify that this initial velocity field has zero vorticity in the inertial frame. We determine the angular 
velocity from Eq. [l], using an iterative method, since velocity derivatives have a weak velocity dependence which 
makes it difficult to find a self-consistent solution directly (see Paper 1, Eqs. A17-A19 and Eq. A23). We find that 
three iterations yields convergence to within a fraction of a percent, and six to within one part in 10 . 

C. Summary of calculations 

We performed several large-scale SPH calculations of NS binary coalescence for the equal-mass case, and a series of 



shorter runs designed to study the effect of varying the mass ratio q. Table III summarizes the relevant parameters of 
all runs performed. Following the notation of Paper 1, we refer to all runs that included both 1PN and 2.5PN effects 
as "PN runs", while those runs performed without 1PN corrections are referred to as "N runs." Note that both N 
and PN runs included the 2.5PN gravitational radiation reaction effects. We did not perform any new completely 



Newtonian calculations (except for a brief test run mentioned in Sec. Ill A) 
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All runs use 5 x 10 4 SPH particles per NS (i.e., the total number of particles N = 10 5 ), independent of q. The number 
of SPH particle neighbors Nn = 100. Shock heating, which is normally treated via the SPH artificial viscosity, was 
ignored, since it plays a negligible role in binary coalescence, especially for fluids with a very stiff EOS (see Paper 1). 
All Poisson equations were solved on grids of size 256 3 , including the space for zero-padding, which yields the proper 
boundary conditions. All PN runs were done using l/c 2 PN = 0.05 and C2.zpn = 2.5, as in Paper 1. Since we have an 
ambiguity in defining the relative time coordinate for different runs, we redefine the initial time to of each run in such 
a way that the time of the first gravity wave luminosity peak tmL = 20, as in run Bl (which starts from t = 0). As a 
result, some runs started before t = (the starting time to for each run is given in the third column of Table pTJ| ). 

Run Al is identical to the N run described in Paper 1, for a binary containing two identical NS modeled as T — 3 
polytropes. This run uses Newtonian gravity with radiation reaction (2.5PN) corrections, and started from an initial 
separation ro = 3.1 R. Runs A2, A3, A4, and A5 feature the same initial conditions, but have mass ratios q — 0.95, 
0.90, 0.85, and 0.80, respectively. Unlike Al, in which the simulation ran until a stable triaxial remnant was formed, 
runs A2-A5 were terminated after the completion of the first peak in the gravity wave luminosity. 

Run Bl is the PN run of Paper 1, also for two identical r = 3 polytropes, but with 1PN and 2.5PN corrections 
included, and an initial separation of ?'o = 4.0-R. Runs B2 and B3, like their Newtonian counterparts, correspond to 
mass ratios q = 0.9 and 0.8, respectively, and the same shorter integration time. 

Run CI is a Newtonian run (with radiation reaction effects) for two identical r = 2 polytropes. Since the ISCO for 
this softer EOS is located at smaller binary separation than for r = 3, (see, e.g., p3[) we start this run with ro = 2.9 R. 
As a test, we compared the results of this run to one with all identical parameters but with a slightly larger initial 
separation, ro = 3.1 R, and found no measurable differences. Run C2 has a mass ratio of q = 0.8 and the same EOS 
and initial separation. Runs Dl and D2 are the 1PN counterparts of CI and C2, with T = 2, an initial separation 
ro = 4.0-R, and q — 1 and q = 0.8, respectively. Run C2 was terminated after the first gravity wave luminosity peak, 
whereas run D2 was continued until a stable remnant configuration was reached, to allow for comparison with Dl. 

All the previous runs started from a synchronized initial binary configuration, as in Paper 1. In contrast, run El 
started from an irrotational initial condition, with the EOS, mass ratio, and initial separation as in run Bl, and 
including all 1PN and 2.5PN terms. It was continued until a stable remnant configuration was reached. 



III. RESULTS 



As in Paper 1, we begin with a brief qualitative discussion of the coalescence process, but focus here on the softer 
EOS with r = 2. In Figs. |l| and 0, we compare the evolution of the system in runs CI and Dl, respectively. Both 
runs are for two identical NS with a T = 2 EOS, but run Dl includes 1PN effects whereas run CI does not. These 
plots can be directly compared to Figs. 1 and 2 of Paper 1, which show the evolution of runs Al and Bl, respectively 
(for two identical NS with T = 3). As seen in Fig. ||, a significant tidal lag angle 8i ag develops in the PN system just 
prior to final merging. Values of 9i ag (calculated as the angle between t he b inary axis and the principal axis of each 
component) upon first contact are given in the fourth column of Table III for all runs. Note that for binaries with 
mass ratios q ^ 1, the secondary develops a larger lag angle than the primary. This effect is much more pronounced 
with the addition of 1PN corrections. In all cases, the merger is accompanied by mass shedding and the formation of 
spiral arms, which may have a significant width. These spiral arms dissipate quickly (by t ~ 40 in Figs. [I] and ^) and 
merge to form a halo of material around the axisymmetric, oblate remnant at the center of the system. It should be 
noted that when we speak of "mass shedding" in this work (both Paper 1 and here), we mean that matter is ejected 
from the central dense core into an outer halo, but not to infinity. By the end of our simulations, all the matter 
forming the outer halo is still gravitationally bound to the system. 



A. Dependence on the NS EOS 

As already pointed out in Paper 1, our PN calculations of NS binary coalescence are most relevant for stiff NS EOS, 
for which most recent calculations give values of GM/Rc 2 ~ 0.1 — 0.2 (for M ~ 1.5 M@). Even if the true NS EOS were 
much softer, making strong GR effects dominant throughout the final binary coalescence, performing hydrodynamic 
calculations in the PN limit would still remain important, since the PN results provide a crucial benchmark against 
which future full-GR calculations can be tested. 

Gravitational radiation wave forms and luminosities for runs Al, Bl, CI, Dl are shown in Figs. || and || The two 
polarizations of gravity waves are calculated for an observer at a distance d along the rotation axis of the system, in 
the quadrupole approximation (see Paper 1, Eqs. 23-25). 
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As was found for NS with a T = 3 EOS (Paper 1), for NS with a T — 2 EOS PN corrections serve to lower the 
maximum gravity wave luminosity. However, regardless of whether we include 1PN terms, the peak gravity wave 
luminosity is larger for the softer EOS. 

Comparing the N and PN runs in Fig. ^, we see that, besides differences in the maximum luminosity, there are 
qualitative differences in the shape of the first luminosity peak. Both PN runs show asymmetry around the peak, with 
a long rise from t ~ 10 — 20 containing a flatter section at t ~ 15 — 20, and then a smooth decrease from t ~ 20 — 25. 
This effect can be attributed to a period of slower binary infall as the two stellar cores first come into contact (as seen 
in Fig. |^, which shows the binary separation r as a function of time). We see that the effect is more pronounced for 
the PN run with aT = 2 EOS (Dl), which has a significant drop in the inspiral rate from t ~ 14 — 18, corresponding 
exactly with a similar plateau in the gravity wave luminosity. The effect is weaker for the PN run with the T = 3 
EOS (Bl), which displays a more symmetric gravity wave luminosity peak. 

Similar behavior is seen in the N runs. For the run with a T = 3 EOS (Al), we see a decrease in the inspiral rate 
from t = 20 — 25, which occurs after the period of maximum gravity wave luminosity, and a simultaneous plateau 
in the luminosity immediately after the peak. For the T = 2 case (CI), we see no apparent decrease in the inspiral 
rate, either before or after t = 20. Correspondingly, this was the only run under discussion which showed an almost 
completely symmetric luminosity peak with respect to time. 

In Fig. y, we show the energy loss to gravitational radiation, calculated as the integral of gravity wave luminosity 
over time. We see that the Newtonian runs, Al and CI, have higher energy losses, in line with the respective peak 
gravity wave luminosities. The secondary peaks in the r = 3 PN run (Bl), discussed extensively in Paper 1, give a 
higher net energy loss than in the T = 2 PN run (Dl), where more energy is lost during the first peak (by a significant 
fraction). Note that energy losses are measured from the initial point of each run. Thus, runs which had negative 
time offsets have E GW > at t = 0. 

More surprisingly, perhaps, we note that a second luminosity peak is visible not only in the T — 2 PN run (Dl), 
but also in run CI, for a T = 2 Newtonian coalescence. In Paper 1, it was shown that the gravity wave luminosity is 
extremely well correlated with the ratio of the first and second principal moments of inertia, I\ and I2 , in effect giving 
a measure of the ellipticity of the remnant in the orbital plane. A similar analysis is shown in Fig. [?], which compares 
runs CI and Dl. Note that we have used a higher density cut for this plot than was used for the T = 3 remnants 
in Paper 1. Here only, we define the inner remnant to consist of all SPH particles with local densities r* > 0.04. 
We see again a strong correlation between ellipticity and gravity wave luminosity. In run Dl, the small- amplitude 
oscillations with period T ~ 5 are caused by an interaction with the outer material of the remnant. If we lower our 
density cut to include all SPH particles with local densities r* > 0.005 (as in Paper 1), thereby including more of the 
tenuous material outside of the core of the remnant, the moment of inertia ratio shows only this oscillation, which 
damps out over time. Since this material contributes only weakly to the quadrupole moment, we conclude that it is 
the dynamics inside the core that controls the gravity wave signal. 

Comparing the results of our run CI, which includes the effects of radiation reaction, with the completely Newtonian 
T = 2 run shown in RS2, we conclude that the asymmetry induced by the larger tidal lag angle is responsible for 
the existence of a second gravity wave luminosity peak, even for systems that eventually reach oblate, non-radiating 
configurations. To confirm this, we checked our results against a completely Newtonian r = 2 test run (without 
radiation reaction terms, and starting from tq = 2.7), and found that a second peak is indeed absent in this case (in 
agreement with RS2). 

Careful inspection of Fig. 7 of Paper 1 (where run Al is referred to as the N run), reveals very small irregularities 
in the otherwise smooth gravity wave luminosity at t ~ 33 and 48, which correspond to the secondary gravity wave 
luminosity peaks found in run Bl. Clearly, the existence of secondary luminosity peaks seems to be reasonably 
universal in these simulations. For run Al, however, modulation of the moment of inertia ratio is virtually absent, so 
that the effect is extremely small in this case. 

B. Dependence on the binary mass ratio 

The dependence of the peak amplitude ft. max = (h\. + h 2 x )max of gravitational waves on the mass ratio q appears 
to be very strong. In RS2, an approximate power law h maK cx q 2 was derived for nearly equal-mass systems, on the 
basis of two purely Newtonian calculations for q — 1 and q — 0.85. This is considerably steeper than the naive scaling 
obtained for two point masses in a Keplerian orbit, which gives /i max oc q. Such a linear scaling is obeyed (only 
approximately, because of finite-size effects) by the wave amplitudes of the various systems prior to final coalescence. 
For determining the maximum amplitude during the merger, however, hydrodynamics must be taken into account. In 
a system with q^l, the more massive star tends to play a far less active role in the hydrodynamics and, as a result, 
there is a rapid suppression of the radiation efficiency as q departs even slightly from unity. For the peak luminosity 
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of gravitational radiation RS found approximately L max oc q 6 . Again, this is a much steeper dependence than one 
would expect based on a simple point-mass estimate, which gives L oc q 2 ( 1 + q) . The results of RS were all for initially 
synchronized binaries, but very similar results have been obtained by Zhuge et al. |pd[ for binaries containing initially 
nonspinning stars with unequal masses. 

The role of the primary and secondary is shown qualitatively in Figs. || and |^ for run D2, which uses aT — 2 EOS 
with q = 0.8, and included 1PN terms. The panels on the left show the primary, center panels the secondary, and those 
on the right the combined system. We see that the evolution is markedly different from that of the q = 1 binary shown 
in Fig. |2j. Here, there is a single spiral arm, which is formed from the secondary as it gets tidally disrupted. Shortly 
after first contact is made, a stream of matter flows from the secondary toward the primary, landing on the trailing 
side (orbital rotation is counterclockwise), as a result of the orbital motion and the significant tidal lags present in 
the system upon contact. As the coalescence proceeds, the secondary is tidally stretched, with the outer portion spun 
out of the system while the inner part is accreted by the primary. Throughout this evolution, the primary remains 
relatively undisturbed, except for a small layer near its surface. By t ~ 40, mass shedding is triggered and the system 
develops a single spiral arm, composed entirely of matter from the secondary as it is completely disrupted. Note that 
this single spiral arm is even wider than those seen in the q — 1 coalescence (with a width comparable to the initial 
NS radius). There is also an extremely small amount of mass shedding from the outer edge of the primary, where it 
joins with a high-velocity stream of matter from the secondary, on the side opposite the single spiral arm. Finally, 
by t ~ 50, what was once the core of the secondary has fallen onto the primary, and the spiral arm has begun to 
dissipate, forming a low-mass halo around the system. 

In Figs. and 0, we show the gravitational radiation wave forms and luminosity for all but two of the synchronized 
binary simulations. For clarity, we only show results for runs Al, A3, and A5 in the r = 3 N plot, since the other 
runs can be safely interpolated from those present. We note that for the r = 3 EOS, the morphology of the peaks 
seen in Fig. ^| seems to be present for all mass ratios. We find plateaus in the gravity wave luminosity before the peak 
for PN runs, and after the peak for N runs. This behavior is much harder to see in the T = 2 runs with q=0.8, which 
have lower peak luminosities in general than their T = 3 counterparts. 



In Fig. 12 we show the maximum gravity wave strain and luminosity for T — 3 and T = 2 EOS binaries, plotted as 
a function of the mass ratio. We find in all cases that the power-law dependence is steeper than would be predicted 
by the point-mass approximation. For the gravity wave strain, we see a slightly steeper power law for N runs than 
for PN runs. For q = 1, N runs have a higher peak strain, but, for both EOS, the q = 0.8 binaries show a higher 
peak strain in the PN case than in the N case. We conclude that the strictly Newtonian scaling obtained by RS2, 
/i mal oc q 2 , remains approximately valid for Newtonian binaries with radiation reaction effects, but is almost certainly 
steeper than the dependence that will be present in realistic NS binaries (based on our results including 1PN effects). 

For the peak gravity wave luminosity, we see a much stronger dependence on the EOS. For the softer, T = 2 EOS, 
we see that the strictly Newtonian RS2 fit of L max cx q 6 (obtained for T — 3 polytropes!) is not nearly as steep as the 
correct relation. Additionally, the difference between N and PN runs in magnitude is less pronounced than for the 
r = 3 EOS, where we have a ~ 20% difference rather independent of the mass ratio. For V — 3, we see that the RS2 
fit is slightly too steep, although the correct power law remains steeper than the point-mass approximation. 



C. Irrotational initial condition 



To study the effects of the initial sp ins o n the evolution of the merger we compare our run El, with the irrotational 
initial condition described in Section II B, with run Bl, the PN run described extensively in Paper 1, which started 
from a synchronized initial condition. Particle plots for run El are shown in Fig. 13, Comparing to Fig. 2 of Paper 1, 
we see that in both cases, the stars develop a large tidal lag immediately before merger, leading to an "off-center" 
collision and a highly asymmetrical merger until t ~ 50. A major difference, however, is the absence of mass shedding 
as the irrotational NS first make contact before t = 20. We see a very small amount of mass shedding and thin 
spiral arm formation from t ~ 25 — 35, involving a very small fraction of the mass that was deposited into the 
outer halo for the synchronized case. Note that mixing of the fluids originating from the two different stars does 
not occur immediately after contact. Instead, as the two NS surfaces come into contact with opposite tangential 
velocities (in the corotating frame), a long vortex sheet forms at the interface between the two fluids (around t ~ 20). 
This vortex sheet is Kelvin-Helmholtz unstable and breaks into a turbulent layer that can then gradually mix the 
fluids pq ] . Unfortunately, the spurious shear viscosity inherent in any discretized simulation is much larger than the 
true viscosity of the NS fluids, and therefore the evolution in the turbulent region is very likely to be dominated by 
resolution-dependent numerical effects. 

Noting that the orbital rotation for the binary is counter-clockwise, we see that all mass shedding into the spiral 
arms actually occurs from the leading edge of each star. Particles that were initially located at small radius relative 
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to the center of mass of the system travel along the vortex sheet between the two stars before being spun off. This 
is in direct contrast to the synchronized case, in which mass shedding occurs almost entirely from regions that are 
initially at large radius, on the outside of each star. Material there is ejected from the merging binary system as it 
spins up and overflows its effective potential well through the outer Lagrange points. Here we find no mass shedding 
via this mechanism, in agreement with previous work on irrotational systems |17| . 

In Figs. [l4| and |l5| we show density contours during the merger of run El, with the velocity field of the fluid in 
the inertial frame overlaid. Mass shedding and spiral arm formation starts at t ~ 22. In Figs. |l6| and |l7| we show 
the density and velocity field of the inner regions of the merger, with the velocities shown with respect to the mean 
corotating frame of the entire system. Specifically, we define the mean corotating angular velocity fi c by taking a 
weighted average over the tangential velocities of all SPH particles, 

fi = J2i m i[( xv y ~~ yvx)/r cy i\i 

where the cylindrical radius is defined as r cy \ = (x 2 + y 2 ) 1 ^ 2 - It is easy to see that this definition reduces to the correct 
orbital angular velocity of the initial binary system, and to the correct angular velocity of any single, rigidly rotating 
object. We see that the vortex sheet formed during the merger lasts from the point of first contact at t ~ 15 until 
t ~ 22. As the coalescence proceeds, turbulence develops along the sheet, with vortices appearing immediately at the 
very center of the remnant, then at several points further out. We note that the central vortex is initially spinning in 
the direction of the orbital motion, as the initial misalignment causes matter to fall in around the center. The outer 
vortices are concentrated around regions of higher density, when matter from deeper inside each NS continues to fall 
onto the vortex sheet. From t ~ 22 to t ~ 25, a slightly different pattern emerges. As the core regions of the original 
NS continue to approach each other, the velocity field between the cores remains turbulent, but a coherent velocity 
field develops around the cores. Counterspinning material located toward the outer parts of the remnant connects 
and forms a simple pattern of circulation. Here the angular velocity appears to become negative as we move away 
from the center (in the inertial frame, this represents the angular velocity decreasing with increasing distance from 
the rotation axis). Material in the center, still somewhat turbulent but retaining the original corotating pattern, has 
much less angu lar momentum, leading to a pattern of strongly differential rotation throughout the merger remnant 



(See Sec. HID) 



The evolution of the mean corotating angular velocity is shown in Fig. |18|. It rises monotonically during the initial 
coalescence until t ~ 25, when the inner NS cores merge. Later it fluctuates with a small amplitude around a nearly 
constant maximum value Q c ~ 0.6 ~ 1.4 x 10 4 rads _1 (M/1.5MQ) 1 / 2 (i?/10km)~ 3 / 2 as the remnant contracts and 
relaxes. By t ~ 50 it settles down with a slowly and steadily decreasing trend as the final triaxial configuration 
continues to lose angular momentum to gravitational radiation. 

The gravitational radiation wave form for run El is shown in Fig. [H| There are no major qualitative differences 
with the wave form obtained for run Bl, which has the same EOS and mass ratio, but a synchronized initial condition 
(see Fig. 3). The gravity wave luminosity is shown in Fig. |2(], along with that of run Bl. We see that irrotational 
initial conditions lead to a slight reduction in the maximum gravity wave luminosity, simply because the system has 
a smaller total angular momentum. The small-amplitude oscillations in the luminosity of run El for t ~ — 15 result 
from small deviations from quasi-equilibrium in the initial binary of magnitude 5p c / p c ~ 0.05. This small amplitude 
gives us confidence that the use of self-similar PN ellipsoids in the construction of our irrotational initial condition 
represents a very small source of error for these calculations. Our irrotational ellipsoids, although still approximate, 
are considerably better than the spherical models used in all previous studies (see, e.g., |TT| , p7| ). 

Both the El and Bl runs show a significant second luminosity peak in the gravity wave emission, although the 
magnitude of the second peak is again smaller for the irrotational case. This second peak occurs at t ~ 33 for both 
runs, indicating that the period between peaks in the gravity wave signal seems to depend primarily on the NS EOS 



and not on the initial spins. Note also in Table [II that the second peaks for runs with aT = 2 EOS are at different 
times, whereas for T = 3, the peaks are always coincident. A third peak is also visible for the irrotational run, but it 
is of much lower amplitude, indicating that the final remnant is more nearly axisymmetric. 



D. Structure of the final merger remnants 

It was found in Paper 1 that the addition of 1PN corrections to the hydrodynamics has the effect of reducing mass 
shedding for mergers of NS with al = 3 EOS. In Fig. |2l|, we show the evolution of the remnant mass M r for runs CI 
and Dl, and find the same behavior for a T = 2 EOS. We show two different density cuts to highlight this behavior. 
By taking all SPH particles with a local density r* > 0.005, we extend our definition of the "remnant" far into the 
outer halo. Instead, a density cut of r* > 0.04 includes just the inner part of the remnant, and excludes any material 
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that was ejected into spiral arms. In both runs mass shedding starts occuring approximately at the time of maximum 
gravity wave emission (t = 20), and lasts for a total time 6t ~ 10, but in the PN case the mass shedding rate is 
significantly smaller, leading to approximately half the total amount of mass shedding as in the Newtonian case. 

In Fig. |22|, we show the radial mass profiles of the inner remnants for the four synchronized q = 1 runs. We see 
that the mass profile at small radii is primarily determined by the EOS, with the r = 2 models showing slightly more 
central concentration, as would be expected, although the inclusion of 1PN effects does decrease the enclosed mass 
in a given cylinder. At r cy \ ~ 1.2, however, we start to see significant differences between runs with and without 
1PN corrections. In the Newtonian runs, matter is ejected much more efficiently to large radii, and thus the final 
mass of the inner remnant falls in the range M r = 1.7 — 1.8, whereas for the PN runs it is M r ~ 1.9. The remaining 
mass, ejected through the spiral arms, forms the halo around the inner remnant. Here we can see from the slopes of 
the mass profiles near r cy \ ~ 2 that for a softer (r = 2) EOS, the density in the halo near the inner remnant is still 
significant, indicating more mass in the inner region of the halo. This is not a surprise since, as already noted by RS2, 
a softer EOS produces wider spiral arms, which in turn dissipate more quickly and transport material less efficiently 
out to very large radii. 

In Table III, we list in the last three columns the total rest mass M r , gravitating mass M gr , and Kerr parameter 
a r of the inner remnant at t = 60. Here the inner remnant includes all SPH particles out to a radius r = 2 (The 
results are rather insensitive to the precise choice of density cut, since the material in the halo is very tenuous). As 
was done in Paper 1, the Kerr parameter for all PN runs is taken to be a r = cJ r /Mg r , where J r is the 1PN angular 
momentum of the remnant. 

A comparison of the angular velocity profiles of the remnants is shown in Fig. ^3l We see, quite surprisingly, that 
the r = 3 run without 1PN corrections (Al) leads to a different form than each of the other runs. All are differentially 
rotating, but in all other cases the angular velocity drops as a function of distance away from the rotation axis, whereas 
in run Al it increases monotonically out to r cy \ ~ 1.4, which is close to the surface of the inner remnant. Only run 
Al agrees with previous purely Newtonian calculations (RSI, RS2), which found that the angular velocity increases 
with increasing distance from the rotation axis. However, in agreement with all previous studies, we find that the 
rotation profiles of all our merger remnants are pseudo-barotopic, i.e., Q is a constant on cylinders. Recent fully 
GR calculations have indicated that differential rotation can increase very significantly the maximum stable mass of 
neutron stars which makes it more likely that merger remnants can be dynamically stable against collapse to a 
black hole. However, differentially rotating configurations could still be secularly unstable on a viscous timescale. 

Run Al is the only one that did not produce a second gravity wave luminosity peak, suggesting that perhaps the 
combination of Newtonian gravity and the stiffer EOS can lead to a lower central angular momentum. In turn, this 
could suppress the quadrupole oscillations shown in Fig. [?]. We also note a general difference in rotational velocity 
between the PN runs and the Newtonian runs, with both N runs showing an increase in angular velocity near the 
surface of the inner remnant, from r cy i ~ 1.2 — 1.4, whereas the PN runs show a steady decrease through this range 
and out into the halo. 

The final structure of the merger remnant appears remarkably independent of the initial NS spins. In Fig. we 
show the mass and angular velocity profiles of the merger remnant for the irrotational binary of run El, compared 
with that of the synchronized system in run Bl. We see that the two merger remnants have very similar profiles. 
Both remnants are differentially rotating with an angular velocity that decreases with increasing distance from the 
axis. The irrotational binary merger remnant has a slightly smaller angular velocity at large radii and a slightly less 
centrally concentrated mass profile. Note also that at large radii, the mass of the remnant for the irrotational case is 
essentially M r — 2, i.e., very little mass shedding occurred. 



IV. SUMMARY AND DIRECTIONS FOR FUTURE WORK 



Using the Lagrangian SPH code described in Paper 1, which is complete to 1PN order and also includes all 2.5PN 
radiation reaction effects, we have further investigated the properties of NS binary mergers in PN gravity. In addition 
to the method presented in Paper 1 for constructing PN equilibrium initial conditions for synchronized binaries, we 
have used the results of the analytic work by LRS based on compressible ellipsoids to create realistic irrotational 
initial conditions. Our approach takes an exact, spherical PN equilibrium model computed for an isolated NS, and 
deforms it according to the shape of the irrotational Darwin-Riemann ellipsoid solution determined from the energy 
variational method of LRS. 

Using the hybrid formalism of Paper 1, in which radiation reaction is treated realistically but 1PN effects are scaled 
down in amplitude to remain numerically tractable, we have studied the effect of PN corrections on several aspects 
of binary NS mergers. We find that NS with a softer EOS, modeled here as T = 2 polytropes, produce higher gravity 
wave luminosities during the merger than NS with stiffer EOS (modeled here and in Paper 1 as T = 3 polytropes). 
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As was found in previous, purely Newtonian calculations, the final merger remnant produced when T — 2 is an oblate 
spheroid, rather than a triaxial ellipsoid (obtained for r = 3). For calculations with and without 1PN effects, we find 
a strong second peak of gravity wave luminosity for a T = 2 EOS, in contrast to the result obtained in Paper 1, where 
we found that for a T = 3 EOS, a strong second peak of emission is present only when 1PN effects are included. 

By holding the EOS and the primary NS mass fixed while varying the secondary NS mass, we have studied 
the dependence of the gravity wave emission on the binary mass ratio. We find that the steep power-law scaling 
h max cx q 2 for the maximum amplitude, derived from purely Newtonian calculations in RS2 for a T = 2 EOS, remains 
approximately correct in PN gravity, for both T = 2 and r = 3 EOS. In both cases the dependence on q appears 
slightly steeper when 1PN effects are neglected (but while retaining 2.5PN effects). The maximum gravity wave 
luminosity for PN mergers follows a slightly steeper power law than the L max oc q 6 found by RS2 for the r = 2 EOS, 
but is significantly flatter for the T = 3 EOS. 

Our coalescence calculation using an irrotational initial condition shows clearly the development of a vortex sheet 
along the surface of contact. A turbulent region is seen until the time when the inner cores of the two NS begin to 
merge, at which point a stable, differentially rotating configuration is created. While the main gravity wave luminosity 
peak is smaller than for a synchronized initial condition, the second peak is of similar relative amplitude and coincident 
in time. A strong third peak was not observed for the irrotational case. The density and rotation profiles of the inner 
remnant appear remarkably independent of the degree of synchronization of the initial binary. 

We find that, in general, the addition of 1PN effects decreases the mass in the outer halo of the merger remnant, 
especially when an irrotational initial condition is used. For this latter case, almost no mass at all is ejected through 
spiral arms during the merger. Since the 1PN corrections are artificially reduced in our calculations, we expect that 
for two NS with realistic parameters and an irrotational initial configuration, no matter at all will be ejected. When 
a synchronized initial condition is used, we find that the halo around the remnant is much denser and less extended 
for a softer NS EOS. 

Our results have potentially important consequences for theoretical models of GRBs based on NS binary mergers 
PUPH - Although current models for long bursts and their associated observed afterglows favor an origin in massive 
star collapse, NS mergers remain an attractive source at least for the separate class of short bursts [|49| . Currently 
the most popular models all assume that a binary NS merger leads to the formation of a rapidly rotating Kerr BH 
surrounded by a torus of ejected material. Energy can then be extracted either from the rotation of the BH or from 
the material in the torus so that, with sufficient beaming, the gamma-ray fluxes observed from even the most distant 
GRBs can be explained [Q. However, our results suggest that the merger of two 1.4 M Q NS with stiff EOS forms an 
object that will not collapse to a BH on a dynamical time (see Paper I, Sec. HID ), in agreement with the preliminary 
full GR calculations of Shibata and Uryu j26j. Moreover, for irrotational binaries (containing two slowly spinning NS; 
see Sec. HI Cj ) , we predict that n o mat ter will be ejected during the merger, implying that no outer torus will form 
around the central core (see Sec. HID ). Even if the central core were to collapse to form a Kerr BH (which is not 
prevented by rotation, since a r < 1; see Table III), energy could not be extracted to power a GRB. 

Our future work using PN SPH will focus on the more realistic irrotational initial conditions introduced in this 
paper, while again varying the other parameters of the problem such as the binary mass ratio and the NS EOS. We 
will also study the dependence of our numerical results on the spatial resolution of the calculations by systematically 
varying the number of SPH particles from N ~ 10 4 all the way to N ~ 10 6 (which will take full advantage of our 
highly optimized parallel SPH code). We will try to establish that our conclusions about the gravity wave sign al and 
final merger remnant configuration are free from numerical artifacts, in spite of the difficulties described in Sec. Ill C . 
Increased resolution should also yield a better understanding of the nature of the vortex sheet formed upon first 
contact of the binary components, and of the details of the transition to a differentially rotating configuration. We 
also plan to study the properties of the gravitational wave signals in Fourier space (as was done for the first time 
by Zhuge et al. |l(]Jll| based on their Newtonian SPH calculations), since the spectral information is crucial to the 
interpretation of future results from LIGO and similar instruments. We will study systematically the dependence of 
the power spectrum of the gravitational radiation on the binary mass ratio, and on the NS EOS, as well as investigate 
the effects of 1PN corrections on the power spectrum. 
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FIG. 1. Evolution of the system for Run CI with T — 2 and q = 1. This run includes Newtonian gravity and radiation 
reaction (2.5PN) corrections. Projections of a random subset of 20% of all SPH particles onto the orbital (x-y) plane are shown 
at various times. The orbital motion is counter-clockwise. Units are such that G = M = R — 1, where M and R are the mass 
and radius of a single, spherical NS. Time has been shifted so that the maximum gravity wave luminosity occurs at t = 20. 
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FIG. 3. Gravity wave signatures for the runs Al (N, F = 3), Bl (PN, F = 3), CI (N, T = 2), and Dl (PN, F = 2). All 
have q = 1. The solid line shows the h+ polarization, the dashed line the h x polarization, both calculated for an observer at a 
distance d along the rotation axis. Note that at t > 60, there is essentially no gravitational radiation given off by T = 2 EOS 
binaries. 
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r=2, q=1.0 
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FIG. 7. Ratio of the principal moments of inertia in the equatorial plane for the remnants of runs CI (N, r = 2; dashed 
lines) and Dl (PN, r = 2; solid lines), compared to the corresponding gravity wave luminosities at late times. We see a 
clear correlation between the two quantities, as the remnants relax towards an axisymmetric, oblate configuration. Here the 
remnants are defined by the density cut r* > 0.04, which includes the entire binary initially, but only the inner part of the 
merger at later times. 
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D2: PN f T=2A q=G,8 




FIG. 8. Evolution of the system for run D2 with F = 2, q = 0.8, and 1PN corrections included. Projections of a random 
subset of 20% of all SPH particles onto the orbital (x-y) plane are shown at various times. Panels on the left show the primary, 
center panels show the secondary, and panels on the right the entire system. 
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D2: PN, T=2.0 f q=0.8 
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FIG. 9. Evolution of run D2, continued. 
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FIG. 10. Gravity wave amplitude h+ for an observer located at a distance d along the rotation axis, comparing systems with 
different mass ratios. The solid lines correspond to q = 1, the dashed lines to q = 0.9, and the dot-dashed lines to q — 0.8. 
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FIG. 12. Dependence of the maximum gravity wave amplitude and luminosity on the mass ratio (for all the synchronized 
binaries). The Keplerian point mass approximation gives ftmax oc q and Z/ max oc q 2 (l + q). Instead, the strictly Newtonian 
hydrodynamic calculations of RS2 give approximate power laws /i max oc q 2 and L m ax oc q 6 for nearly equal-mass binaries 
containing T = 3 polytropes. Note that both axes are plotted logarithmically. 
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El: PN, T=3 J O i q=1.0 J Irrotational 
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FIG. 13. Evolution of the V = 3, q = 1 irrotational binary system (El). Conventions are as in Fig. |l|. We see evidence for 
spiral arms forming after the initial merger, but containing much less mass than in the synchronized case. Note also that a 
longer time is required after merger to relax to a triaxial configuration. 
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FIG. 14. Density contours and velocity field for the irrotational system in run El. Density contours are spaced logarithmically, 
4 per decade. The velocity field is shown in the inertial frame. A vector representing the speed of light (c = C2.5PN, the physical 
value) is shown for comparison. 
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FIG. 15. Density contours and velocity field for run El, continued. 
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El: FN, r=S.O, q=LO» Irrotational 



0.5 - 



-0.5 

-1 
1 

0.5 



-0.5 - 



-1 




-1 -0.5 0.5 1 - 



FIG. 16. Density contours and velocity field for the irrotational system in run El, here focusing on the inner region of the 
merger. Panels correspond to the same times as shown in Fig. |l4|. Density contours are spaced logarithmically, 8 per decade. 
The velocity field is shown in the corotating frame , as defined by Eq. [l|. As the surfaces of the two stars come into contact, 
they form an vortex sheet which extends over the full length of the interface. 
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El: FN, r=2.0, q=1.0, Irrotational 
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FIG. 17. Density contours and velocity field in the inner region for run El, continued. Panels correspond to the same times 
as in Fig. 15. 
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El: PN, r=3.0, q=1.0, Irrotational 
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FIG. 18. Evolution ol the mean angular velocity of the corotating frame fi c for run El (as defined by Eq. [|). 
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El: PN, T=3, (i 1.0. Irrotational 




31 



PN, r=3, q=1.0 




20 40 60 

t 



FIG. 20. Gravity wave luminosity for run El (irrotational) compared with that of run Bl (with a synchronized initial 
condition). Both runs include 1PN corrections, and have F = 3 and q = 1. 
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T=2, q=1.0 




FIG. 21. Evolution ol the remnant mass for the T = 2, q = 1 runs (CI and Dl). Suppression of mass shedding in the PN 
case leads to a higher final remnant mass. The density cut r* > 0.04, shown as dashed lines, includes the inner remnant only. 
The cut r, > 0.005, shown as solid curves, extends further into the outer halo. 
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FIG. 22. Radial mass profiles of the final merger remnants in runs Al, Bl, CI and Dl (at t = 60). Here r cy \ is the distance 
from the rotation axis. We see that the internal structure is governed primarily by the EOS, with the Y — 2 models slightly 
more centrally condensed than the r = 3 models, although 1PN effects do decrease the enclosed mass at small radii. At larger 
radii, we find more mass contained in the PN remnants, since much less mass has been ejected through spiral arms. 
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FIG. 23. Angular velocity profiles for the same merger remnants shown in Fig. 22. We see that all are differentially rotating, 
and all but the Newtonian T = 3 remnant show a decrease in angular velocity with increasing radius for r cy \ < 1.3. In general, 
at larger radii, the PN models show slower rotation, regardless of the EOS. 
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PN, r=3, q=1.0 
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FIG. 24. Comparison of the merger remnants from runs Bl (synchronized initial condition) and El (irrotational initial 
condition) . Both runs have T = 3 and 1PN corrections. Both the mass profiles and the rotation profiles of the inner remnants 
appear to be rather independent of the initial NS spins. 



36 



TABLE I. Properties of the r = 2 polytropic NS sequence. Units are denned such that G = M = R = 1. The polytropic 
constant k is defined by P — fcr*, and fcjy is the Newtonian value, which is used in the Lane-Emden equation. The central 
density (r*) c is in units of M/R 3 , while the central pressure ratio (P/r*c 2 ) c and central gravitational potential ratio (U*/c 2 ) c 
are dimensionless. 
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TABLE II. Properties of the PN polytropic models for NS of varying masses. 
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TABLE III. Selected quantities from each of the simulations. Note that only runs Al, Bl, CI, Dl, D2 and El were extended 
until a stable final remnant was formed. All other runs were terminated shortly after the first gravity wave luminosity peak. 
Here to is the time at which the run was started (see Sec. II C), and is the lag angle at first contact, given for both NS (see 
Sec. pi] ). Quantities involving the first and second gravity wave luminosity peaks are labeled with superscripts (1) and (2). All 
SPH particles with r < 2.0 are us ed to calculate the final rest mass M r , gravitational mass M gr , and Kerr parameter a r of the 



final remnant (at t — 60; see Sec. HID) 
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